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Glueballs and Hybrids (Gluons as Constituents) 

Doug Toussaint ^ 

'^Department of Physics, University of Arizona, Tucson, AZ 85721, USA 



After a brief introduction to hybrid and glueball source operators, I summarize recent lattice results for these 
particles. 
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In addition to the familiar baryons made from 
three quarks and mesons made from a quark 
and an antiquark, the spectrum of QCD is ex- 
pected to contain glueballs, hybrids and multi- 
quark states. Since the earliest lattice QCD cal- 
culations it has been recognized that lattice QCD 
should be used to calculate the properties of these 
particles. Hopefully, lattice results will help to 
sort out which of the many known hadrons con- 
tain important glueball or hybrid components. 
The experimental situation for glueballs and hy- 
brids with light quarks was recently reviewed in 
Ref. jl]. Evidence for an exotic 1 ^ particle has 
been seen in several experiments!^, and is now 
substantial enough to merit inclusion in the Re- 
view of Particle Properties Unfortunately, the 
mass of this particle, 1400 MeV, is quite a bit 
lower than the lattice results so far, and its decays 
conflict with phenomenological flux-tube models. 

Of course, we would like to know everything 
about these particles — their masses, their de- 
cays, their wavefunctions, their mixings .... 
However, the place to start, and where most of 
the calculations have stopped, is with a calcula- 
tion of the masses. This works just like spectrum 
calculations for ordinary hadrons: find an oper- 
ator with the desired quantum numbers Oi{x,t), 
and compute 
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The ground state is the easiest to find; for ex- 
cited states you must manipulate several Oi and 
look for combinations with A^^^ = 0, or simulta- 
neously fit to several masses. Even for the ground 



state, it is still important to construct an operator 
minimizing overlap with excited states. 

To make these operators, we can combine the 
following ingredients: 

q"--. quark, color 3, 

q"". antiquark, color 3, 

q'^Tq'^: color singlet quark bilinear 

q°-Tq^: color octet quark bilinear 

B"^: color magnetic field, color 8, J^^ = 1+" 

E°'^: color electric field, color 8, J^^ = 1^^ 

L: orbital angular momentum 
For example, to make the simplest glueballs you 
might use: 

BfBf = 0++ 

BfBf - ^Trace{) = 2++ 



One way to understand the construction of the 
simplest hybrid operators is to replace one of the 
glue fields in these glueball operators with a qq 
bilinear. In particular, take a quark and in the 
IS state, "break" the color, and insert the color 
magnetic field. In its simplest form, called "Ba- 
sic" in the upper left of Fig. ^ the color magnetic 
field is just the difference between the product 
of links in a small loop traversed clockwise and 
counterclockwise connecting the quark spinor to 
the antiquark. 

The very first thing you might try is to take a 
pion operator, 9°759°, and combine it with the 
color magnetic field to get q'^'j^q^Bf''. The quan- 
tum numbers of this operator are found by com- 
bining the quantum numbers of the pion and the 
magnetic field — ^ 1^ =1^^ — the same 
quantum numbers as the p. 

Similarly, if you take the quark and antiquark 
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to have spin one, q'^jq"', you can combine it with 
the color magnetic field to produce three different 
quantum numbers: 1 = 0"+®l~+©2~+ 

Q-+. qa^.qbj^ab. num. of pion) 

e^Jkq°'liq^Bf■. (exotic ) 

2-+ 



q°'^iq Bj +{i ^ j) — Tr(): (quantum num- 
bers of TT2) 

Thus, if hybrid spectroscopy works like con- 
ventional meson spectroscopy or glueball spec- 
troscopy, we expect that the lowest multiplet of 
hybrid states contains 1^^, 1 ^ and 2 ^ 
particles, both isovector and isoscalar, which will 
be split by the color hyperfine interactions. The 
most interesting of these is the 1 ^, which has 
exotic quantum numbers. 

If we continue, and consider qq operators with 
the quark and antiquark in a P wave in the nonrel- 
ativistic limit, we can construct operators includ- 
ing exotics with quantum numbers 0^ and 2"' , 
and by using the color electric field you can make 
a — exotic {ai E = 1++ (g) 1 — = — . . .). 

The "Basic" cartoon of a hybrid operator in 
Fig. 1^ is not really acceptable since it violates the 
cubic symmetries of the lattice — there is nothing 
special about the upper left plaquette. An obvi- 
ous way to symmetrize it is to average over all the 
nearby plaquettes, producing the "clover" in the 
upper right of Fig. |l|. Another way to proceed is 
to begin by moving the traversal of the plaquette 
in one of the two directions down to the lower left, 
as in the center of Fig. |l|. This doesn't give the 
full cubic symmetry, but (with a few more loops 
out of the plane of the figure) does have rotational 
symmetry around one axis. Then one can imag- 
ine pulling the quark and antiquark apart along 
this axis, reaching the cartoon in the lower right 
of Fig. 1^ by way of the intermediate stage in the 
lower left. In this way we make contact with a 
picture of a hybrid meson as an excited flux tube 
connecting the quark and antiquark. 

For both glueballs and hybrids it is necessary 
to use "smeared" or "fuzzed" links for the gluon 
parts of the operator to get good overlap with 
the ground state particles. The details of how 
this is done vary, and more often than not the 
parameters are determined empirically. 



Basic 



rot. (cubic) sym. 



rot. sjrtn. around 1 axis 



pull quarks apart 





Figure 1. Construction of hybrid operators, il- 
lustrating the connection between the "local" op- 
erators and "excited potential" operators. 



Along with the choice of the operator which 
creates the hybrid, we can choose how the quark 
and antiquark propagate. The ideal way is to use 
the full relativistic formulation. However, this is 
time consuming, and for heavy enough quarks it 
is much more sensible to use the nonrelativistic 
(NRQCD) propagators. In fact, one can go all 
the way to the infinite mass limit, and simply 
use static quarks. (This was the first method 
used to study hybrids on the lattice). In ad- 
dition to its simplicity, the use of static quarks 
gives a very physical picture of the hybrid me- 
son as a quark and antiquark moving in a poten- 
tial that depends on the state of the fiux tube 
connecting them (keywords: adiabatic potential, 
Born-Oppenheimer approximation), and within 
this approximation it is straightforward to com- 
pute the wavefunctions of the hybrids, as well as 
masses of excited states. (Since the quarks can't 
move in the static approximation, it is necessary 
to use the "flux tube" versions of the operators 
in Fig. 1^ so that the quark and antiquark can be 
created at the desired separation.) 

While the first lattice calculation of hybrids was 
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done in 1983, progress was slow until the last 
few years. Highlights of these calculations can 
be found in Refs. through |]l9| . 

All of the lattice studies agree with our expec- 
tation that the 1 ^ should be the lightest exotic 
hybrid, and this particle has received the most at- 
tention. To give a flavor for the state of the art. 
Table |^ shows a selection of results for the mass 
of the 1 hybrid or, in the case of heavy quarks, 
for the mass difference between the 1 ^ hybrid 
and the qq S-wave mass. The tabulated results 
show a pleasing convergence with time, and for 
recent results (bold face) the agreement among 
different groups is quite satisfactory. 

Glueball mass calculations have a longer his- 
tory than hybrid mass calculations. Most have 
been done in the quenched approximation, per- 
haps partly because glueballs are the actual ex- 
citations of quarkless QCD. The state of the art 
is represented by three accurate calculations in 
the quenched approximation: a UKQCD calcula- 
tion[Q, a series of works by the GFll group 
and a recent calculation by Morningstar and 
Peardon[^. The first two of these are conven- 
tional calculations, using isotropic lattices with 
the one-plaquette action, while the third uses 
anisotropic lattices with an improved action in 
the spatial planes. All three use results at sev- 
eral lattice spacings to get good control of the 
a ^ limit for the lowest glueballs. In addition, 
the anisotropic calculation [|2| gets convincing re- 
sults for a few of the excited states. As empha- 
sized by Teper[||, the UKQCD and GFll cal- 
culations of the 0++ glueball mass in units of the 
lattice spacing are in excellent agreement, and the 
difference in their reported results for the mass in 
MeV comes from different methods of assigning a 
lattice spacing. Some of this ambiguity is inher- 
ent in quenched calculations, since we expect that 
quenched mass ratios in the a —^ limit will differ 
from real world mass ratios. To compare these 
calculations I follow contemporary practice and 
use To to set the lattice spacing. Ref. quotes 
masses in units of ro, and for the conventional 
action calculations I convert to physical units us- 
ing the interpolating formula for vq of Guagnelli, 
Sommer and Witting using rp — 0.50 fm. 
With this definition of lattice spacing, some of 



the results of these three calculations are: 
UKQCD 

0++ 1645(50) MeV 

2++ 2337(100) MeV 

GFll 

0++ 1686(24)(10) MeV 

2++ 2380(67)(14) MeV 
Morningstar and Peardon 

0++ 1659(43)(16) MeV 

2++ 2304(8) (24) MeV 

0-+ 2494(28)(24) MeV 

0++* 2561(173)(28) MeV 

2++* 3499(43) (35) MeV 
All three calculations are in reasonable agree- 
ment for the O"^"*" and 2++ glueball masses, and 
the Morningstar and Peardon results show the 
state of the art in extracting excited state masses. 

There is also a growing body of work, much of it 
by the GFll group [2l|], on the decay rates of glue- 
balls, and on the mixing of glueballs with quarko- 
nia to produce the physically observed hadrons. 
Here the idea is to take the quenched glueball 
masses and the quenched light quark isoscalar 
and ss quarkonia masses, together with either 
a measurement or a model for their overlap, 
and diagonalize a matrix to find the physical 
states. (Oddly, this is a calculation where we 
actually prefer a quenched mass spectrum to a 
full QCD lattice calculation.) The eigenstates are 
then identified with the experimental /o(1370), 
/o(1500) and /,/(1710) (assuming J = 0). Such 
an analysis has been done by the GFll group us- 
ing their own numerical results pT[| and by Teper 
using numerical inputs from the literature p3|. 
Teper proposes: 

/(1710) = 0.42g + 0.90ss + 0.13wM 
/(1500) = 0.77.g - 0.43ss + 0.48mm 
/(1370) -0.49.9 + 0.10SS -I- 0.87?2m 

while the GFll group suggests: 

/(1710) = 0.86(5)5 + 0.30(5)ss-h0.41(9)uM 
/(1500) = -0.13(5).g + 0.91(4)ss-0.40(ll)?2M 
/(1370) = -0.50(12)5 + 0.29(9)ss -I- 0.82(9)mm 

where g is the glueball component, ss the 
strangeonium component, and uu the (isoscalar) 
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light quark component. 

These results are quite different. The largest 
difference conies from differences in the unmixed 
masses, with Teper having rriuu = 1.36, rriss = 
1.61 and rrig = 1.48 GeV and the GFll group 
having rriuu — 1A7, rngs = 1.51 and irig = 1.63 
GeV. Note that in the first analysis the unmixed 
ss is heavier than the glueball, while in the sec- 
ond it is lighter. Not surprisingly, the heavi- 
est physical particle is mostly ss or mostly glue- 
ball, respectively. Clearly it is important to 
know whether the quenched ss is heavier than 
the quenched glueball, and the GFll group has 
attacked this question by calculating quarkonium 
and glueball masses on the same samples. 

It is also interesting to test the effects of dy- 
namical quarks on the glueball spectrum. In prin- 
ciple this is tricky, since a full QCD spectrum cal- 
culation with a glueball source operator will pro- 
duce the masses of the physical states, which are 
mixtures of glueballs and quarkonia. In practice, 
the calculations that have been done so far have 
used quark masses large enough that the quarko- 
nium mass exceeds the glueball mass, and it is 
reasonable to simply say thay they have measured 
the glueball mass. The largest calculation is from 
the TxL/SESAM collaborations [|5|, and another 
preliminary result from UKQCD was presented at 
The TxL/SESAM results for 
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the 0^"*" glueball arc within errors of the quenched 
results, although they do see hints of a larger de- 
pendence on the lattice size. However, the pre- 
liminary UKQCD results have a much smaller 
mass, despite being done at approximately the 
same sea quark mass, as measured by rriT^/mp. 
The UKQCD results are done on much coarser 
lattices, albeit with the clover action rather than 
the Wilson quark action. (UKQCD used (3 = 5.2 
while TxL/SESAM used 5.6.) The 0++ glueball 
mass has long been known to be small on coarse 
lattices|^, but the preliminary UKQCD results 
are even smaller than we would expect from our 
experience with the quenched theory. 

Since the exotic 1 ^ signal found in experi- 
ments at 1400 MeV is much lower than expected 
from lattice calculations (and most other theo- 
retical approaches), it is tempting to ask whether 
it could be something else, most likely a 4-quark 



Table 1 

Some results for 1 

Date Ref. 
hhg - bb: 
1990 il 
1993 [10| 
1997 pi 



1997 [15| 

1997 [14 1 

1998 [18| 

1999 jl9| 
ccg — cc: 
1990 M 
1996 [12| 

1998 [17| 

1999 |l8| 
ssg 

1996 Inl 
1996 12 1 
udg 

1996 Inl 

1996 [12| 

1998 [16| 

1998 121 
Abbreviations: 
= Relativistic, An. 
Clo. = clover. 

Notes: (1): value with a determined differently, 
(2): a = 0.095 fm, (3): a = 0.075 fm, (4): Model 
to extrapolate to niq — 0, 120 MeV below ss 
mass, a = 0.095 fm, (5): Extrapolation from 
several aniq values, a — 0.075 fm, (6): Nf = 2 
dynamical quarks, extrapolate from several amq 
values, a — 0.086 fm, (7): Same as (5). 

(qqqq) state. In principle, this question is an- 
swerable with lattice methods, but it is a difficult 
subject. Nonetheless, there is a small but grow- 
ing body of work on lattice 4-quark states [^^ , 
beginning with the simplest case where all four 
quarks are static, and moving into the case of two 
static and two moving quarks. An amusing limit 
has recently been studied by Michael and Penna- 
nen|29|, where the two quarks are very heavy, and 
the two antiquarks light. The two quarks have an 
attractive interaction in the 3 color combination, 
and since they are very heavy they can bind into 



Method 


AM (GeV) 




St. 


1.11(3)(?) 




NR. 


0.8(?)(?) 




NR. 


1.68(10) 






1.40(14) 


(1) 


NR. 


1.14(21) 




St. 


1.3 




NR. (An.) 


1.542(8) 




St+NR(An.) 


1.49(2)(5) 




St. 


0.94(3) 




Rel.(Wil.) 


1.34(8)(20) 




Rel.(Clo.) 


1.22(15)(?) 




NR. (An.) 


1.323(13) 






M (GeV.) 




Rel.(Clo.) 


2.00(20) 


(2) 


Rel.(Wil.) 


2.17(8)(20) 


(3) 


Rel.(Clo.) 


1.88(20) 


(4) 


Rel.(Wil.) 


1.97(9)(30) 


(5) 


Rel.(Wil.) 


1.90(20) 


(6) 


Rel.(Clo.) 


2.11(10)(?) 


(7) 



St. = static, NR. = NRQCD, Rcl. 
anisotropic, Wil. = Wilson, 
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a small object which is effectively an antiquark. 
Then the three antiquarks may be regarded as an 
antibaryon, with one constituent very heavy. 

I expect that lattice calculations will continue 
to be important in sorting out the rich structure 
of the QCD spectrum. It is important to extend 
our work to calculations of decay rates and mix- 
ings among particles. 
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